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ABSTRACT 

Context. MGRO J2019+37 is an unidentified extended source of very high energy gamma-rays originally reported by the Milagro 
Collaboration as the brightest TeV source in the Cygnus region. Its extended emission could be powered by either a single or several 
sources. The GeV pulsar AGL J2020.5+3653, discovered by AGILE and associated with PSR J2021+3651, could contribute to the 
emission from MGRO J2019+37. 

Aims. Our aim is to identify radio and near-infrared sources in the field of the extended TeV source MGRO J2019+37, and study 
potential counterparts to explain its emission. 

Methods. We surveyed a region of about 6 square degrees with the Giant Metrewave Radio Telescope (GMRT) at the frequency 
610 MHz. We also observed the central square degree of this survey in the near-infrared A" s -band using the 3.5 m telescope in 
Calar Alto. Archival X-ray observations of some specific fields are included. VLBI observations of an interesting radio source were 
performed. We explored possible scenarios to produce the multi-TeV emission from MGRO J20 19+37 and studied which of the 
sources could be the main particle accelerator. 

Results. We present a catalogue of 362 radio sources detected with the GMRT in the field of MGRO J2019+37, and the results 
of a cross-correlation of this catalog with one obtained at near-infrared wavelengths, which contains ~ 3 x 10 5 sources, as well as 
with available X-ray observations of the region. Some peculiar sources inside the ~1° uncertainty region of the TeV emission from 
MGRO J2019+37 are discussed in detail, including the pulsar PSR J2021+3651 and its pulsar wind nebula PWN G75.2+0.1, two 
new radio-jet sources, the Hn region Sh 2-104 containing two star clusters, and the radio source NVSS J202032+363158. We also 
find that the hadronic scenario is the most likely in case of a single accelerator, and discuss the possible contribution from the sources 
mentioned above. 

Conclusions. Although the radio and GeV pulsar PSR J2021+3651 / AGL J2020. 5+3653 and its associated pulsar wind nebula 
PWN G75.2+0. 1 can contribute to the emission from MGRO J2019+37, extrapolation of the GeV spectrum does not explain the 
detected multi-TeV flux. Other sources discussed here could contribute to the emission of the Milagro source. 

Key words, gamma-rays: observations - H n regions - infrared: stars - radio continuum: stars - X-rays: binaries 



1. Introduction 

The Galactic very-high-energy (VHE) y-ray sources discovered 
by the latest generation of Cherenkov observatories (H.E.S.S., 
MAGIC, Milagro) are currently an actively studied topic in mod- 
ern high-energy astrophysics. Among the ~75 detected sources, 
nearly one third remain yet as unidentified. A significant num- 
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ber of them have extended morphologies on 0.1-1° scales in the 
TeV energy band, ensuring that the identification of counterparts 
at lower energies is a very difficult task. The most representative 
of this new population of Galactic sources is TeV J2032+4130, 
inside whose error box both compact and extended radio sources 
on arcsecond scales were found (Paredes et al. 2007). XMM- 
Newton observations of this source also detected faint extended 
X-ray emission (Horns et al. 2007a). 
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An addition to the population of extended, unidentified TeV 
sources was reported by the Milagro collaboration, following the 
discovery in the Cygnus region of the most extended TeV source 
known so far (Abdo et al. 2007 a,b). The TeV emission from this 
area covers several square degrees and includes diffuse emission 
and at least one new source, MGRO J2019+37, located to within 
an accuracy of ±0.4°. After the Crab Nebula, MGRO J2019+37 
is the strongest source detected by Milagro. The Tibet AS-y ex- 
periment confirmed the detection of this source by measuring a 
5.8cr signal compatible with the position of MGRO J20 19+37 
(Amenomori et al. 2008). On the other hand, VERITAS inferred 
an upper limit that is compatible with the Milagro detection for 
a hard-spectrum extended source (Kieda et al. 2008). 

The origin of all these types of emission and their associa- 
tion with astrophysical sources is unclear. Although a possible 
connection with the anisotropy of Galactic cosmic rays was pro- 
posed (Amenomori et al. 2006), the TeV y-ray flux measured at 
12 TeV from the diffuse emission of the Cygnus region (after ex- 
cluding MGRO J2019+37), exceeds that predicted by a conven- 
tional model of cosmic ray production and propagation (Abdo et 
al. 2007a). This strongly infers the existence of hard-spectrum 
cosmic-ray sources and/or other types of TeV y-ray sources in 
the region. It is unclear whether the emission originates in ei- 
ther a single extended source or a combination of several point 
sources. MGRO J2019+37 is positionally coincident with the 
EGRET sources 3EG J2021+3716 and 3 EG J2016+3657 (see 
Fig. 1). These sources could represent the GeV counterparts 
to the TeV source MGRO J2019+37, which may be a multi- 
ple source. Only one of them, 3EG J2021+3716, appears in the 
bright gamma-ray source list published by the Fermi Gamma- 
ray Space Telescope (Abdo et al. 2009). Previous observations 
with AGILE illustrated its pulsar nature and inferred an associa- 
tion of this source with PSR J202 1+3651 (Halpern et al. 2008). 

To explain steady VHE y-ray emission, hadronic models 
have been developed by several authors (e.g., Aharonian & 
Atoyan 1996; Butt et al. 2003; Torres et al. 2004; Bordas et al. 
2009). The electromagnetic radiation produced by both hadronic 
jets from microquasars and Galactic cosmic rays, and their in- 
teraction with the ISM were explored by Bosch-Ramon et al. 
(2005). The interaction between the high energy protons, accel- 
erated at the jet termination shock, and the interstellar hydrogen 
nuclei produces charged and neutral pions (n~, n + and 7r°); the 
first set will decay to electrons and positrons and the second set 
to photons. The primary radiation, 7r°-decay photons, is in the 
y-ray band, but the secondary particles can produce significant 
fluxes of synchrotron (from radio frequencies to X-rays) and 
bremsstrahlung emission (from soft y-rays to the TeV range), 
and in general lower efficiency, inverse Compton (IC) emis- 
sion by interaction with ambient infrared photons. Detectable 
fluxes of extended and steady emission should be produced by 
this mechanism. Other scenarios involve a jet-driven termina- 
tion shock at which relativistic electrons produce synchrotron 
and TeV IC emission (Aharonian & Atoyan 1998). In this con- 
text, X-ray observations provide a crucial constraint of the IC 
emission. 

To understand the nature of the Milagro source in the Cygnus 
region, we performed a multiwavelength campaign comprising 
a deep radio survey at 610 MHz using the Giant Metrewave 
Radio Telescope (GMRT) interferometer covering the 3.5°x3.5° 
MGRO J2019+37 field, near-infrared observations in the K s 
band using the 3.5 m telescope at Calar Alto of the central square 
degree, and archival X-ray data. 

This paper is organized as follows. In Sect. 2, we report on 
previous radio surveys of the Cygnus region, while in Sect. 3 we 



present the GMRT survey and the results obtained. In Sect. 4, we 
provide an overview of the near-infrared survey, and in Sect. 5 
we report on the cross-correlations both between our GMRT sur- 
vey and the near-infrared survey, and between the GMRT survey 
and previous X-ray observations. We comment on particularly 
interesting sources in Sect. 6 and we discuss their possible con- 
tribution to the TeV emission of MGRO J2019+37 in Sect. 7. 
We finish with our conclusions in Sect. 8. 

2. Previous radio surveys of the Cygnus region 

At radio frequencies, the Cygnus region has been imaged many 
times, sometimes as part of Galactic surveys. However, these 
studies were carried out at poor angular resolution and/or a 
relatively high limiting flux density. Some of the most repre- 
sentative of previous surveys are: the Canadian Galactic Plane 
Survey (CGPS) performed with the Synthesis Telescope at the 
Dominion Radio Astrophysical Observatory (DRAO) at 408 and 
1420 MHz, with angular resolutions of 5.'3 and 1.'6, and limiting 
flux densities of 9 and 1 mJy, respectively, at declination of +40° 
(Taylor et al. 2003); the Westerbork Synthesis Radio Telescope 
(WSRT) 327 MHz survey with an angular resolution of 1' and 
a limiting flux density of 10 mJy (Taylor et al. 1996); and the 
DRAO 408 and 1430 MHz survey with angular resolution of 
3.5' x 5.2' and 1.0' x 1.5', respectively, and limiting flux densities 
of 150 and 45 mJy, respectively (Wendker et al. 1991). The most 
recent survey of this region is the WSRT 350 and 1400 MHz 
continuum survey of the Cygnus OB2 association, with angular 
resolutions of 55" and 13", and limiting flux densities of 10-15 
and 2 mJy, respectively (Setia Gunawan et al. 2003). The WSRT 
survey does not cover the MGRO J2019+37 field. 

3. GMRT 610 MHz Radio Survey 

3.1. Observations 

The MGRO J2019+37 region was observed with wide-field 
deep radio imaging at 610 MHz (49 cm) using the GMRT, lo- 
cated in Pune (India). We designed an hexagonal pattern of 19 
pointings to cover the region of about 2.5°x2.5° centred on the 
MGRO J2019+37 peak of emission. The observations were car- 
ried out in July 2007, but were affected by a series of power 
failures in the array and compensatory time was scheduled in 
August 2007. 

The flux density scale was set using the primary amplitude 
calibrators 3C 286 and 3C 48, which were observed at the be- 
ginning and end of each observing session. On the other hand, 
phase calibration was performed by repeated observations of the 
nearby phase calibrator J2052+365. Each pointing was observed 
for a series of scans to achieve a good coverage in the uv plane, 
the total time spent on each field being 45 minutes. The total 
effective time amounts to 20 hours. 

Observations were made in two 16-MHz upper and lower 
sidebands (USB and LSB) centered on 610 MHz, each divided 
into 128 spectral channels. The data of each side-band were sep- 
arately edited with standard tasks of the Astronomical Image 
Processing System (MPS) package. There were no major radio 
frequency interference (RFI) problems. However, we did find 
that narrow band RFI affected a few channels across the band, 
which were completely flagged. Once poor antennas, baselines, 
or channels were removed, the bandpass correction was used to 
extend the calibration to all channels. After the bandpass calibra- 
tion, the central channels of each sideband were averaged, lead- 
ing to a data file of 5 compressed channels, of a bandwidth small 
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Fig. 1. Radio map obtained with the GMRT at 610 MHz (greyscale) convolved with a circular restoring beam of 30". The red cross and box indicate 
the center of gravity and its positional uncertainty including statistic and systematic errors of the TeV emission from the source MGRO J2019+37 
(Abdo et al. 2007b). The conspicuous radio sources located inside this box correspond to the extended H n region Sh 2-104, a bright compact radio 
source also detected with the VLA as NVSS J202032+363158, and two newly discovered jet-like sources (A and B). The position probability 
contours (50%, 68%, 95%, and 99%, from inside to outside) of the Third EGRET catalogue sources 3EG J2021+3716 and 3EG J2016+3657 
(Hartman et al. 1999), as well as the GeV source GeV J2020+3658 (blue ellipse) (Lamb & Macomb 1997) are superimposed. The magenta cross 
indicates the position of the pulsar wind nebula PWN G75. 2+0.1. The blazar B2013+370 within 3EG J2016+3657 is also labeled. 



enough to avoid bandwidth smearing problems in our images. 
Standard calibration for continuum data was performed beyond 
this point. At the end of the self-calibration deconvolution iter- 
ation scheme, we combined both USB and LSB images of each 
pointing and mosaicked the entire region using the AIPS task 
FLATN. 



We produced different maps of between high and low angu- 
lar resolution of the GMRT mosaic. Our highest quality image 
has an rms of 0.2 mJy beam -1 with a 5" resolution because of 
the long baselines of the GMRT. A low angular resolution ver- 
sion was also produced using a restoring beam of 30" to enhance 
the extended radio sources in the field. This map has an rms of 
0.5 mJy beam -1 . 



3.2. Results 



Figure 1 shows a low angular resolution radio image of the 
MGRO J2019+37 field, together with the position of sources 
at other wavelengths. The location of MGRO J2019+37 is con- 
sistent with those of the EGRET sources 3EG J2016+3657 
and 3EG J2021+3716. The first of them is positionally coin- 
cident with the blazar-like source B2013+370 (G74.87+1.22) 
(Mukherjee et al. 2000; Halpern et al. 2001), although this 
blazar is well outside the inner box of MGRO J2019+37. The 
second is marginally coincident with the pulsar wind nebula 
PWN G75.2+0.1 (Hessels et al. 2004). High-energy gamma-ray 
pulsations originating in the pulsar were detected by AGILE and 
Fermi (Halpern et al. 2008; Abdo et al. 2009). There are other 
known strong and/or extended radio sources in the field, such as 
the brightest one inside the MGRO J2019+37 center of gravity 
box, NVSS J202032+363158, and the Hn region Sh 2-104 (also 
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Fig. 2. Number of sources versus log 5^ eiik for the 362 sources detected 
in the GMRT 610 MHz radio survey. 



known as Sh 104). Other interesting sources not obvious at first 
glance become evident when considering the whole field in de- 
tail. Some of them display a resolved morphology, and in Sect. 6 
we discuss these objects in more detail. 

3.3. Radio catalogue 

We applied the SEXtractor automatic procedure (Bertin & 
Arnouts 1996) to our 5" resolution mosaic (with pixel size of 
1") to produce a list of sources with peak flux density higher 
than about ten times the local noise after primary beam correc- 
tion. Objects with less than 5 connected pixels above the thresh- 
old were not included. The output was visually inspected and 
all candidate detections inferred to be false (i.e., deconvolution 
artifacts near bright sources) were simply deleted by hand. We 
used the local background analysis in SEXtractor to take into ac- 
count the uneven background because of beam response effects. 
Considering the 5" beam size of the mosaic that we used, and the 
signal-to-noise ratio that we required for detection, we estimate 
that the positions obtained have a typical uncertainty of 075 or 
smaller. 

The resulting list, considered to be very reliable although 
not complete at the lowest flux density levels, contains 362 ra- 
dio sources. Among them, 203 are fainter than 10 mjy and the 
majority were previously undetected at radio wavelengths. We 
present the catalogue in Table 2 of the online material accom- 
panying this paper. The first and second columns provide the 
catalogue number and the source name. The third and fourth 
columns give the J2000.0 position in right ascension order. The 
fifth and sixth columns provide the peak flux density and the 
local noise, respectively. The seventh and eighth columns list 
the integrated flux density and its error. Uncertainties quoted for 
the peak and integrated flux densities are based on the formal 
errors of the fit and allow the reliability of the detection to be 
judged. However, they do not include the error due to primary 
beam correction as a function of angular distance to the phase 
centre because of unknown antenna offsets, which is estimated 
to be around 10% of the flux density values (see for instance 
Paredes et al. 2008). In Fig. 2, we show the source distribution 
histogram as a function of log S ^ eak . 



4. Near-infrared survey 

We also carried out a near-infrared (NIR) survey of the cen- 
tral square degree of the region using the OMEGA2000 wide 
field camera (15' x 15') on the 3.5 m telescope at Centra 
Astronomico Hispano Aleman (CAHA) in Calar Alto (Spain) 
on 25 September 2007. This instrument consists of a Rockwell 
HAWAII2 HgCdTe detector with 2048 x 2048 pixels sensi- 
tive from 0.8 to 2.5 /xm. The observations were performed in 
the /T s -band (2.15 yum) to minimize the interstellar absorption. 
Individual frames were sky-subtracted, flat-field corrected, and 
then combined into a final mosaic using the AIPS task FLATN. 
The ensamble of 4 X 4 pointings covers almost completely the 
center of gravity and uncertainty region of the TeV emission 
from the source MGRO J2019+37. The average limiting mag- 
nitude across the mosaic is K s 17 mag, and the total field of 
view is 0.9° x 0.9°. Astrometric solutions for the final frames 
were determined within ±07 1 by identifying about twenty refer- 
ence stars in each pointing, for which positions were retrieved 
from the 2MASS catalogue (Skrutskie et al. 2006). A cata- 
logue of ~315 000 near-infrared sources was produced using the 
SEXtractor package. 



5. Radio, near-infrared, and X-ray cross-correlation 
catalogue 

We performed a cross-correlation between the radio and near- 
infrared source catalogues. Considering the 075 uncertainty in 
the radio positions and the 071 uncertainty in the NIR ones, 
we used a conservative maximum offset of 076 for associa- 
tions (neglecting systematics between both catalogues). There 
are 42 of the 362 detected radio sources inside the area im- 
aged in the near infrared. A total of 6 of these 42 sources have 
a near-infrared counterpart candidate within 076 of their radio 
position. Their magnitudes are listed in the ninth column of 
Table 2 of the online material. The chance coincidence prob- 
ability of finding a NIR source closer than 076 to a given ra- 
dio source is estimated to be the number of NIR sources mul- 
tiplied by the area of the uncertainty in positions occupied by 
the 42 radio sources, divided by the total area of the region: 
(315 000 x 42 n 076 2 )/(3240" x 3240") = 1.4. Therefore, of 
the six radio sources with NIR counterpart, we expect that one 
of them is a random coincidence. 

We also obtained source lists of all X-ray observations of 
the region performed by Chandra and XMM-Newton, computed 
by the celldetect and edetect_chain tasks from CIAO 4.0 
and SAS 8.0, respectively. A total of 41 of the 362 radio sources 
are located in fields observed in X-rays, which cover an area of 
314 arcmin 2 (1 130973 arcsec 2 ) and contain 519 X-ray sources. 
We found that 5 of the 41 radio sources have an X-ray coun- 
terpart candidate within 5" (the typical uncertainty for XMM- 
Newton). Their X-ray fluxes are listed in the tenth column of 
Table 2 of the online material. The chance coincidence prob- 
ability of finding a radio source closer than 5" to a given X- 
ray source is estimated to be the number of radio sources mul- 
tiplied by the area of the uncertainty in positions occupied by 
the 5 X-ray sources, divided by the total area of the region: 
(41 x 519 7T 5" 2 )/(l 130973 arcsec 2 ) = 1.4. Therefore, of the 
five X-ray sources with radio counterpart we also expect that 
one of them is a random coincidence. 

A single triple radio/near-infrared/X-ray coincidence has 
been found (source number 115 in Table 2 of the online ma- 
terial). 
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6. Individual sources in the MGRO J2019+37 field 

The most interesting radio sources that appear in the uncertainty 
region of the TeV emission (red box in Fig. 1) are described 
below. 

6.1. PSR J2021+3651 1 PWN G75.2+0. 1 

The radio pulsar PSR J2021+3651 has a rotation period P = 
0.104 s, a characteristic age of P/2P = 17 kyr, and a spin- 
down luminosity E = 3.4 x 10 36 erg s _1 . It is coincident with 
the unidentified source GeV 2020+3658 (Roberts et al. 2002), 
which overlaps with the EGRET source 3EG J2021+3716 (see 
Fig. 1). Chandra observations of this pulsar detected a ~ 20" x 
10" pulsar wind nebula named PWN G75 .2+0.1 (Hessels et al. 
2004). Chandra observations of the pulsar and its PWN detected 
rings and jets around PSR J202 1+3651, and inferred a distance 
to the pulsar of 3^1 kpc (Van Etten et al. 2008), in contrast to 
the 12 kpc implied by the pulsar dispersion measure (Roberts et 
al. 2002). XMM-Newton observations show emission extending 
to a distance of ~ 10-15 arcminutes, whereas radio observations 
with the VLA at 1 .4 GHz show a radio nebula coincident with 
the X-ray extension (Roberts et al. 2008). 

AGILE detected the source AGL J2020. 5+3653 at energies 
above 100 MeV range, which shows pulsations and was asso- 
ciated with the pulsar PSR J2021+3651 (Halpern et al. 2008). 
The photon spectrum of the source can be fitted by a power-law 
of photon index r = 1.86 ± 0.18 in the range 100-1000 MeV, 
while a turndown is seen above 1 .5 GeV. This source appears as 
1AGL J2021+3652in the first AGILE catalog of high confidence 
gamma-ray sources (Pittori et al. 2009). Fermi also detected the 
source 0FGL 2020.8+3649 in positional coincidence with the 
pulsar (Abdo et al. 2009). 

We found neither a (low-frequency) radio nor a near-infrared 
source at the position of PSR J2021+3651. The nearest near- 
infrared source is at a distance of 3'.'9 and has a K s magnitude 
of 17.3. In the radio, from our 610 MHz GMRT data we can 
establish an upper limit to any possible point-like counterpart 
of 1.0 mJy by multiplying the background emission level by 
a factor of 5. The radio flux density of the extended emission 
found with the VLA at 1 .4 GHz amounts to ~700 mJy in an area 
of about 100 arcmin 2 , which for a uniform distribution yields 
7 mJy arcmin 2 . On the other hand, the rms of our low-resolution 
radio map at 610 MHz shown in Fig. 1, with a beam size of 30", 
is 0.5 mJy beam -1 . This provides a conservative 5-cr upper limit 
of either 2.5 mJy beam -1 or 9 mJy arcmin -2 . This upper limit im- 
plies that if the radio emission is uniformly distributed, its spec- 
tral index must be above -0.3. This value is compatible with the 
radio emission being produced by the synchrotron mechanism, 
as expected in this nebula. 

6.2. Jet-like radio sources 

We discovered two jet-like radio sources located well inside the 
uncertainty region of MGRO J2019+37. Their J2000.0 positions 
are a=20 h 18 m 32 s , 5=+37°02.'5 (source A) and a=20 h 19 m 48 s , 
5=+37°06.'7 (source B). In Fig. 3, we show a GMRT high 
resolution image of each of them superimposed on the near- 
infrared image. Both sources appear to be unresolved in the 
NVSS 1.4 GHz catalogue (Condon et al. 1998). Based on our 
GMRT survey, the NVSS survey, and the Westerbork Northern 
Sky Survey (WENSS; Rengelink et al. 1997) data at 327 MHz, 
we estimate a spectral index of -1.16+0.02 for source A, clearly 
indicating a non-thermal nature. It is interesting to note that the 



source is not detected in the VLA Low-Frequency Sky Survey 
(VLSS; Cohen et al. 2007) at 74 MHz, with a 3-cr upper limit 
of 1.2 Jy. With the spectral index above, we would expect a 
flux density of 1.9 Jy, clearly indicative of a turnover at lower 
radio frequencies, which could be produced by intrinsic self- 
absorption or by Galactic foreground free-free absorption. The 
GMRT and NVSS data for source B provide a spectral index 
of -0.7+0.6, compatible with the non-detection in WENSS, and 
suggesting a non-thermal nature for this radio source. 




2018 35 34 33 32 31 30 29 

RIGHT ASCENSION (J2000) 




2019 50 49 48 47 46 

RIGHT ASCENSION (J2000) 



Fig. 3. Radio and near-infrared image composition of the jet-like ra- 
dio sources A (top) and B (bottom). The GMRT radio contours are 
superimposed on the A" s -band 3.5 m CAHA telescope images. Top: 
Source A. Contours correspond to 5, 9, 15, 23, 45, 60, and 80 times 
0.16 mJy beam -1 , the rms noise. The integrated flux density of the 
source is 164.7+0.3 mJy. Bottom: Source B. Contours correspond to 
5, 8, 12, 20, 30, and 42 times 0.16 mJy beam -1 , which is the rms noise. 
The integrated flux density of the source is 28+0.1 mJy. The synthe- 
sized radio beams of 5" are plotted in the lower-right corners of both 
images. 

Source A (Fig. 3-top) shows a double-sided morphology, 
sources #141 and #142 in Table 2 of the online material, with 
a slight bending towards the south-east. This structure resem- 
bles ones typically seen in radio galaxies with a non-negligible 
pressure from the intergalactic medium. Unfortunately, there is 
no clear extended NIR counterpart in the axis joining the radio 
lobes that could be identified with the parent galaxy, and no firm 
conclusion can be obtained from the present data. 

Source B (Fig. 3-bottom) shows a morphological and spec- 
tral similarity to the radio lobes of the 'great annihilator' 
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Fig. 4. Top: Composite radio and near-infrared image centred on the Sh 2-104 region. The contours correspond to 10, 20, 35, 55, 80, 100, 
125, 155 times the rms noise of 0.3 mjy beam -1 of our GMRT 610 MHz (49 cm wavelength) image. We overlay our A^ s -band near-infrared 
image of the same field obtained using the 3.5 m telescope at Calar Alto. The blue cross marks the position and 1-cr uncertainty of the 
ROSAT source 2RXP J201742. 3+3645 13. Bottom-left: Young massive stellar cluster deeply embedded in a UCHII region found on the east- 
ern rim of Sh 2-104. Bottom-right: New cluster candidate in the center of Sh 2-104 previously assumed to be a single star (identified later as 
2MASS J20174184+3645264). The scale of each image is indicated by the horizontal bar. The colour scale of the bottom images was changed to 
display the individual stars within each cluster more clearly. 
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IE 1740.7-2942, a microquasar at the Galactic center (Mirabel 
et al. 1992). The two lobes correspond to sources #193 and 
#194 in Table 2 of the online material. We did not detect 
a radio core in this source but, as for the one present in 
IE 1740.7-2942, it could have a flat spectrum and the flux 
density at such a low frequency is expected to be very low 
compared to that of the radio lobes. As can be seen in the 
figure, there are two near-infrared objects close to the cen- 
tral position of the source. Their J2000.0 coordinates and mag- 
nitudes are: a=20 h 19 m 47 s :74, £=+37°06'4072, K S =16.5, and 
a=20 h 19 m 47 s :86, <5=+37°06'3979, K S =17A. Their proximity sig- 
nificantly biases the photometry. The bright source is point-like 
and offset from the axis traced by the radio lobes. The faint 
source is aligned with the axis and fuzzy, implying that the origin 
of the double radio source is most likely a radio galaxy. 

Previous ROSAT pointed observations (Obs. Id. 500248P 
conducted on 24 October 1993) did not detect any of these 
two radio-jet sources, placing a 3-<x upper limit of 7 x 10~ 14 
erg cm 2 s _1 on their persistent flux in the energy range 0.1- 
2.4 keV. With the present data, we cannot elucidate whether the 
sources are Galactic or extragalactic, although there are hints of 
their extragalactic nature. 

6.3. Hn region Sh 2-104 

Sh 2-104, also known as Sh 104, is an optically visible Hn re- 
gion of 7' diameter at a distance of 4.0+0.5 kpc (Deharveng et al. 
2003). There is a central 06 V star suspected of being responsi- 
ble for ionizing the region (Lahulla et al. 1985). The appearance 
of Sh 2-104 in the optical and in the radio bands is very similar, 
although the radio images show the presence of an ultra compact 
H ii (UCHII) region at the eastern border, which is not visible in 
the optical image (Deharveng et al. 2003). The interaction be- 
tween the expanding H n region Sh 2-104 and the UCHII region 
may be responsible for triggered star formation in the latter, re- 
sulting in a deeply embedded young cluster. This region has also 
been detected as a high luminosity (3 x 10 4 L ) IRAS source. 

Our GMRT observations (see Fig. 4) detect a structure sim- 
ilar to that found at 1.46 GHz with the VLA (Fich 1993) and at 
1.4 GHz within the NVSS radio continuum survey (Condon et 
al. 1998). 

We also observed Sh 2-104 in the near-infrared A" s -band. The 
images obtained are deeper than those from 2MASS. Figure 4 
shows our near-infrared images of the field of Sh 2-104 with the 
radio emission contours superimposed. In the eastern region of 
the ring (to the left side), the near-infrared image shows the well 
known cluster associated with the UCHII region, which must 
contain at least one massive OB star (Deharveng et al. 2003). In 
the central part of the image, the single 06 V star of Lahulla et 
al. (1985), which corresponds to 2MASS J20174184+3645264 
(Skrutskie et al. 2006), now appears to be resolved as sev- 
eral point-like objects, indicative of the presence of a clus- 
ter. Therefore, apart from this ionizing early-type star, the new 
cluster candidate could also contribute to the formation of 
the Hn region (e.g., additional early type stars, wind shocks). 
Furthermore, an elongated arc along the south of Sh 2-104 as 
well as to the east of the UCHII region can be discerned in the 
NIR images. These features may be related to the interaction be- 
tween the expanding H n region and the interstellar medium. 

Despite its deep coverage at other wavelengths, Sh 2- 
104 was poorly explored in the X-ray domain. Previous X- 
ray observations of this region by ROSAT detected a source 
(2RXP J201742.3+364513; Rosat consortium 2000) located at 
Qr= 20 h 17 m 42 s .3, <5=+36°45'13" with a positional error of -12", 



Table 1. Non-simultaneous flux density measurements of the source 
NVSS J202032+363158 obtained from different surveys. 



Survey 


Frequency 


Flux density 




(MHz) 


(mJy) 


VLSS 


74 


6354+708 


WENSS 


327 


1442+216 


CGPS 


408 


1180+360 


GMRT 


610 


833+ 56 


NVSS 


1400 


386+ 12 


GB6 


4850 


121+ 12 


87GB 


4850 


108+ 15 



overlapping with the central star 2MASS J20 174 184+3645264 
and the cluster candidate (see bottom-right of Fig. 4). The count 
rate of (4.1 ± 0.5) x 10~ 3 count s _1 in the energy range 0.1- 
2.0 keV, provides a flux of (5.8 ± 0.7) x 10~ 14 erg cirr 2 s" 1 
based on the assumption of a thermal spectrum with a tem- 
perature of 1.5 keV (a typical value for colliding wind re- 
gions). On the other hand, OB stars are known to be X-ray 
sources, presumably because of shocks in their stellar winds 
(see Giidel 2004 for a review). According to the complete study 
by Berghofer et al. (1997) of more than 200 isolated OB stars 
detected in ROSAT data, for an 06 V star, with bolometric lu- 
minosity of 8 x 10 38 erg s _1 (Martins et al. 2005), the corre- 
sponding X-ray luminosity is 1.2 x 10 32 erg s _1 . Considering 
a distance of 4.0 kpc to both Sh 2-104 and its ionizing central 
star 2MASS J20 174 184+3645264, the expected X-ray flux is 
6.3 x 10~ 14 erg cirr 2 s _1 , fully compatible with the detected X- 
ray flux from 2RXP J201742.3+364513. 

6.4. NVSS J202032+3631 58 

The source NVSS J202032+363 158 is the brightest compact ra- 
dio source within the error box of the TeV peak emission of 
MGRO J2019+37 in our GMRT observations. Its flux density 
at 610 MHz is 833 mJy and has not been resolved. This source 
appears in the VLSS at 74 MHz, in the WENSS at 327 MHz, 
in the CGPS at 408 MHz (Taylor et al. 1996), in the Effelsberg 
survey of the Cygnus X region at 1420 MHz (Wendker et al. 
1991), and in the Green Bank 4.85 GHz northern sky surveys 
87GB (Gregory & Condon 1991) and GB6 (Gregory et al. 1996). 
In Table 1, we summarize the detected flux densities within 
these surveys, and we plot the corresponding spectrum in Fig. 5. 
Assuming a stable flux density, the radio spectrum of this source 
can be described by S v = (523 ± 2) mJy [y/GHz]- a94±001 , and it 
is, therefore, clearly a non-thermal emitter. Despite very low fre- 
quencies being sampled, no evidence of the turnover frequency 
below 1 GHz is obvious from this simple power-law fit. 

By inspecting of the NRAO archives, we found a previous 
VLA snapshot (6 min on source) of this radio source at the 
20 cm wavelength in B configuration (providing a nominal syn- 
thesized beam of 4") observed on 25 March 1989. This observa- 
tion was calibrated using standard AIPS tasks, including phase 
self-calibration. A uniformly weighted image is shown in Fig. 6- 
left. As can be seen, this radio source is resolved, displaying 
a one-sided radio jet extending a few arcsec towards the north, 
with a core component of ~250 mJy and a secondary component 
of about 70 mJy. To enhance the compact structure of the source, 
we obtained an image for the longest baselines of the same VLA 
run, using a Mv-range of 30-50 kA. The image, shown in Fig. 6- 
right, clearly shows a compact core and a secondary component, 
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Fig. 6. (Left): Image of the source NVSS J202032+363158 at 21 cm obtained using uniform weights on B-configuration VLA data. The source 
is resolved, displaying a one-sided radio jet extending a few arc-sec towards the north, with a core component of ~250 mjy and a secondary 
component of about 70 mjy. The rms of the image is 0.22 mjy beam -1 . Contours correspond to 4, 8, 16, 32, 64, 128, 256, 512, and 1024 times 
the rms noise. (Right): Image from the same data performed using an wv-range of 30-50 k/l, which clearly shows the core and the component 
discussed in the text, with peak flux densities of 170 and 20 mjy, respectively. The rms of the image is 2.5 mjy beanT 1 . Contours correspond to 
-3, 3, 4, 6, 8, 14, 20, 30, 40, 50, and 60 times the rms noise. The synthesized radio beams are plotted in the lower-right corners of both images. 
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Fig.5. Radio spectrum of NVSS J202032+363158 based on the flux 
densities compiled in Table 1 . The straight line is a simple power-law 
fit. 

with peak flux densities of 170 and 20 mjy, respectively, resem- 
bling the large-scale jet of a microquasar. 

To explore the source at higher angular resolutions, we ob- 
served the core of NVSS J202032+363158 at 1.6 GHz (18 cm 
wavelength) with the European VLBI Network in eVLBI mode 
(eEVN). This is a technique in which the signals from distant 
radio telescopes are directly streamed into the central data pro- 
cessor for real-time correlation, instead of being recorded on 
disk or tape. The observation took place on 3 March 2007 from 
5:00 to 13:00 UT (centered on MJD 54163.375), and was per- 



formed using 6 antennas: Cm, Mc, Jb-2, On-85, Tr, and Wb. 
Scans on NVSS J202032+363158 were interleaved with scans 
on the compact phase calibrator J2015+3710, with a 6-min cy- 
cle time (66 s on the calibrator and 246 s on the source). The 
data were recorded using dual polarization and 2-bit sampling, 
at 256 Mbps. A total bandwidth of 32 MHz per polarization was 
provided by 4 sub-bands. The e- VLBI data were processed at the 
Joint Institute for VLBI in Europe (JIVE) correlator in real time, 
using an integration time of 2 sec. The target source was cor- 
related with the position obtained from the VLA-B (30-50 k/l) 
data: a J20 oo.o = 20 h 20 m 33 s .0401 and <5 J20 oo.o = +36°31'57'.'480, 
for a total maximum uncertainty of 100 mas. During observa- 
tions, we experienced synchronization problems and the corre- 
lation had to be restarted several times. A few antennas were 
dropped out of the correlation jobs during the gaps used for mea- 
suring the system temperatures. Due to these disconnections, 
part of the data, which is not recorded onto disks for these ex- 
periments, was lost during the correlation, and the true on-source 
time is estimated to be around 3 hours. 

We performed the post-correlation data reduction using the 
MPS software package and Difmap. We applied ionospheric 
corrections to the visibility data, and the system temperatures 
were used to obtain the a priori visibility amplitude calibration. 
All stations produced fringes with the 1-Jy phase calibrator, sit- 
uated at 1°2 from the target, and we therefore transferred the so- 
lutions for the phases to the target source. We improved the am- 
plitude calibration using correction factors for each antenna ob- 
tained from the self-calibration of J2015+3710. Self-calibration 
of the NVSS J202032+363158 data was impossible because of 
the lack of bright sources in the primary beam of the antennas. 
The phased-referenced natural-weighted image that we obtained 
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had a synthesized beam of 22.7x19.8 mas at a position angle of 
30°3, and an rms noise of 0.20 mJy beam -1 . No significant de- 
tections were found within a distance of 5" from the correlated 
phase center. 

There is no near-infrared counterpart candidate to 
NVSS J202032+363158. The nearest sources are both at 
4'.'1, with magnitudes of 14.2 and 17.3 in the Zf s -band. 

7. Could any of the selected individual sources 
power the TeV emission from MGRO J2019+37? 

MGRO J2019+37 covers a sky region of approximately 1° x 1°. 
The extended emission could be produced by either a single 
powerful accelerator, or by the superposition of several point- 
like sources. Although we focus on the individual sources pre- 
sented in Sect. 6, we cannot exclude some of the additional radio 
sources listed in Table 2 of the online material being responsible 
for, or contributing to, the Milagro source. 

If MGRO J2019+37 is a single extended source, and not 
a combination of different sources, the origin of the >12 TeV 
emission is likely to be hadronic. The time required to fill a re- 
gion of a size of ~1° (or (l-5)xl0 20 cm at 2-10 kpc distance) 
with electrons of ~ 100 TeV by means of diffusion is 

fdiiF = 1.5xl0 12 ^ B_ 6 s , (1) 

where /?2o = R/10 20 cm is the source size, and = B/IQ' 6 G 
is the ISM turbulent magnetic field. For realistic ISM densi- 
ties of «ism < 10 4 cm -3 , and reasonable magnetic/mm-far IR 
field energy densities, i.e., >1 eV cm 4 , the electron cooling 
timescale is dominated by synchrotron and IC losses and found 
to be f coo i < 10 11 s. Therefore, given that fdiff » fcooi> electrons 
injected from a single accelerator cannot fill the entire multi-TeV 
source. Otherwise, protons cool mainly by means of collisions 
with the ISM nuclei (pp) 

fcool ~ 10 15 /«ISM S. (2) 

For n ISM ~1 000-200 cm -3 (2-10 kpc) in the Milagro region, 
proton injection luminosities of ~ 10 37 erg s -1 should be enough 
to explain the observed luminosities (Abdo et al. 2007a) as- 
suming that ~0.1-1% of the proton power is in >12 TeV pho- 
tons. As a result of pp interactions, secondary electron-positron 
pairs and neutrinos are also produced with luminosities and en- 
ergies similar to those of gamma-rays (e.g., Kelner et al. 2006). 
These secondary pairs should radiate via synchrotron, relativistic 
bremsstrahlung, and IC. Extended radio and X-ray emission was 
detected within the Milagro region (Roberts et al. 2008; Hessels 
et al. 2004; Van Etten et al. 2008). However, the smaller extent 
of these diffuse sources compared to the size of the TeV emis- 
sion makes any possible association difficult. Nevertheless, for 
typical ISM densities and magnetic and radiation fields, most of 
the emission from the secondary pairs could be produced at rel- 
atively low gamma-ray fluxes, rendering them undetectable. 

Once the most probable emission scenario is decided, we 
will be able to see whether the different objects proposed in 
Sect. 6 could act as the accelerator. 

The spin-down luminosity of PSR J202 1 +365 1 is marginally 
in agreement with the energetic requirements stated above. 
Nevertheless, for this object to act as the accelerator, most of this 
luminosity should be in the form of protons (as in, e.g., Horns 
et al. 2007b). In addition, the accelerated protons should escape 
the ~ 10' nebula in a time shorter than or equal to the age of 
the pulsar, w 17 kyr, which may not be possible if the turbu- 
lent magnetic field in the nebula reaches value of several 10 pG 



or higher. On the other hand, the turndown in the AGILE GeV 
spectrum questions the association of AGL J2020.5+3653 as the 
only counterpart to MGRO J2019+37. Extrapolation of the last 
two data points in the spectrum shown in Halpern et al. (2008) 
provide a flux at 20 TeV a factor of 3500 below the reported 
MGRO J2019+37 flux (Abdo et al. 2007b). Even ignoring the 
turndown and fitting the entire spectrum with a single power-law, 
there is still a one order of magnitude difference. Therefore, if no 
additional components are present in the GeV- TeV spectrum of 
PSR J2021+3651/PWN G75.2+0.1, this source alone can hardly 
explain the multi-TeV emission from MGRO J2019+37. 

The massive star and the star-forming region (MSR; SFR) 
associated with the Hn region Sh 2-104 could be responsible 
for the extended Milagro source if they were capable of in- 
jecting ~10 37 erg s -1 in the form of relativistic protons into 
their surroundings. Assuming an efficiency of a 10% for the 
kinetic energy converted to non-thermal proton energy in the 
shocks present inside the Hn/SFR region, about 100 mas- 
sive (proto)stars producing jets or winds with velocities of ~ 
10 8 cm s" 1 and mass-loss rates of ~ 10~ 6 M yr~' would be 
required to reach the needed proton luminosities. It seems un- 
likely that Sh 2-104 can harbor such a high number of massive 
(proto)stars. However, Sh 2-104 may be part of a larger MSR 
or SFR that have not yet been detected, and in that case, the 
larger whole MSR or SFR could represent the emitter of the 
whole Milagro source through wind collisions or jet/medium in- 
teractions (see, e.g., Torres et al. 2004; Romero 2008), respec- 
tively. In this scenario, thermal free-free radio emission from the 
whole SFR would be expected. The non-detection of this SFR in 
our GMRT observations could be explained by free-free absorp- 
tion in the ionized regions and the surrounding material of the 
MSR/SFR, although the development of the particular details of 
this scenario are beyond the scope of this work. Observations 
searching for maser emission with instruments such as Apex or 
Nanten could help us to detect this hypothetical star-forming re- 
gion. 

We note that the accelerator itself might be outside 
MGRO J2019+37, as in the case of the stellar cluster 
Berkeley 87 mentioned in Abdo et al. (2007a). This cluster could 
accelerate the protons that would then escape from it diffusing 
towards, and ultimately interacting with, a denser region located 
near the Milagro source best-fit model position. 

We found three non-thermal radio sources with jet-like struc- 
tures in the field of MGRO J2019+37: sources A and B, and 
source NVSS J202032+363158. Although some arguments sup- 
port the extragalactic nature of sources A and B, we cannot ex- 
clude the possibility that they are Galactic in nature. VLBI ob- 
servations of NVSS J202032+363 158 provide an upper limit of 
1 mJy beam -1 to the flux for a beam size of ~20 mas. Therefore, 
this source did not exhibit a compact core during our obser- 
vations. It was either completely resolved or is a variable ra- 
dio source, since no radio emission is expected in the high/soft 
state of microquasars. In any case, these three sources could be 
hadronic microquasars whose jets would interact with the ISM 
accelerating protons (e.g., Heinz & Sunyaev 2002; Bordas et al. 
2009). The accelerated protons may escape from the accelerating 
region colliding with the surrounding regions of the ISM, render- 
ing very high-energy emission (e.g., Bosch-Ramon et al. 2005). 
From the energetic point of view, although these sources could 
explain the Milagro source, the lack of clear X-ray counterparts 
needs to be explained, if accretion is taking place in these ob- 
jects. It might be the case that accretion is inefficient in produc- 
ing X-rays (as could be the case in LS 5039; e.g., Bosch-Ramon 
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et al. 2007). Finally, a microquasar located outside the Milagro 
region could be powering the multi-TeV radiation. 

The constraint on a hadronic origin for the Milagro emis- 
sion does not apply if the source consists of different accelera- 
tors/emitters. In this case, several leptonic emitters, which may 
or may not coincide with (some of) the sources discussed here, 
could be behind MGRO J2019+37. 

From this analysis, we conclude that several objects 
should be considered when trying to understand the origin of 
MGRO J2019+37, although the nature of the accelerator/emitter 
remains uncertain. Insights into this question could be provided 
by further multiwavelength studies, by future imaging atmo- 
spheric Cherenkov telescopes (MAGIC-II, H.E.S.S.-II), and by 
Fermi, which should be able to constrain the source position 
and morphology more tightly, and explore in detail its physics 
by obtaining spectral information across a broad wavelength 
range. Finally, neutrino detections with future neutrino instru- 
ments could provide additional evidence to support the hadronic 
scenario. 



8. Conclusions 

We have carried out a deep radio survey of about 6 square de- 
grees region in the direction of MGRO J20 19+37, and a near- 
infrared survey of the central square degree. This has provided a 
catalogue of 362 radio sources and a catalogue of 315 000 NIR 
sources. The radio and NIR data presented here detect a large 
number of previously unknown sources and shed additional light 
on known objects. We have found that if a single accelerator is 
powering MGRO J2019+37, the most likely origin of the multi- 
TeV emission is hadronic in nature. We have shown that the 
extrapolation of the spectrum of the pulsar AGL J2020.5+3653 
does not explain the detected flux from MGRO J2019+37. This 
indicates either that there is an additional component in the GeV- 
TeV spectrum of the pulsar and/or that other sources, such as 
those discussed here, could contribute to the emission of the 
Milagro source. The results presented in this paper may be use- 
ful in interpreting future data provided by the Fermi satellite of 
the gamma-ray sources in this remarkable region of the Galactic 
plane. The physical understanding of the most relevant sources 
in the field is currently a work in progress, in addition to the 
analysis of new XMM-Newton and AGILE observations. 
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Table 2. List of GMRT sources detected at 610 MHz, including their names, positions, peak flux densities, integrated flux densities, NIR mag- 
nitudes (dots denote that the radio source is outside our NIR mosaic, while ND represents non-detections) and X-ray fluxes (dots imply that the 
radio source is outside the X-ray fields, while ND stands for non-detections). 
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j.17-18-78 6 
t J / . Jo.ZO.O 
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0.42 


7.03 
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ND 
ND 
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20:27:12.44 


+36:58:18.8 
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GMRTJ202733.3+373116 
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GMRT J202735.0+373124 


20:27:35.00 


+37:31:24.6 


5.64 



0.26 25.51 0.10 

0.43 10.25 0.12 

0.32 12.67 0.15 

0.16 10.10 0.10 

0.34 7.58 0.14 

0.27 10.33 0.13 ... ND 

0.25 5.82 0.11 ... 1.9 ± 0.4 x 10~ 14 

0.25 5.82 0.09 

0.29 5.94 0.09 

0.20 2.52 0.06 

0.24 6.07 0.08 

0.30 41.84 0.20 

0.32 57.57 0.14 

0.32 9.50 0.16 

0.48 99.53 0.21 

0.48 10.46 0.19 

0.38 89.73 0.20 ... ND 

0.50 5013.70 1.40 ... ND 

0.32 14.10 0.17 

0.30 25.96 0.30 

0.31 13.32 0.18 
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